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Abstract 

The Quadratic Assignment Problem (QAP) is a well-known NP-hard 
combinatorial optimization problem that is at the core of many real-world 
optimization problems. We prove that QAP can be written as the sum 
of three elementary landscapes when the swap neighborhood is used. We 
present a closed formula for each of the three elementary components and 
we compute bounds for the autocorrelation coefficient. 



1 Introduction 

We will define a landscape for a combinatorial problem using a triple (X, N, /), 
where / : X (-» R defines the objective function and the neighborhood operator N 
assigns a set of neighboring solutions N(x) S X to each solution x. If y <E N(x) 
then y is a neighbor of x. The landscape that is induced can be used as a search 
space for optimization using local search. 

There is a special kind of landscape which is of particular interest due to 
their properties. They are the elementary landscapes, and are characterized by 
the following equation: 

avg{/(y)} = /(z) + §(/-/(*)) (1) 

y£N(x) a 

where d is the size of the neighborhood, |iV(x)|, which we assume the same 
for all the solutions in the search space, / is the average solution evaluation 
over the entire search space, and A; is a characteristic constant. Equation ([TJ 
is usually called Grover's wave equation and makes it possible to compute the 
average value of the fitness function / evaluated over all of the neighbors of x; 
we denote this average using &vg{f(y)} yeN ( x y 



avg{/(y)} = — -— /(») 

y£N(x) PWI yeN(x) 



Other properties also follow. Assuming f(x) ^ / then 



f{x) < min < avg{/(y)}, / \ V f(x) > max < avg{/(y)}, / 

yeN(x) y€N(x) 



This implies that all maxima are greater than / and all minima are less than / 
[4] . A landscape (X, N, /) is not always elementary, but even in this case it is 
possible to characterize the function / as the sum of elementary landscapes [3] , 
called elementary components of the landscape. 

The Quadratic Assignment Problem (QAP) is an NP-hard combinatorial 
optimization problem 2 . A lot of research has been devoted to analyze and solve 
the QAP. Some other problems can be formulated as special cases of the QAP. 
One important example is the Traveling Salesman Problem (TSP). The QAP 
is not an elementary landscape when the swap neighborhood is considered pQ. 
But, to the best of our knowledge the exact expressions for the elementary 
components of the QAP are not known. 

Such decomposition could be useful from the theoretical and practical points 
of view. In theory, the landscape decomposition of QAP can be used to compute 
the exact expression of the autocorrelation functions and the autocorrelation co- 
efficient pp. In practice, the landscape decomposition together with the Grover's 
wave equation can be used to compute the average value of the objective func- 
tion in the neighborhood, which can be used as a base for new operators or 
algorithms. In particular, a new family of selection operators can be designed 
which select the individuals according to the average fitness value in the neigh- 
borhood of a solution x instead of using the fitness value of the solution itself. 
These selection operators could be especially useful to distinguish solutions that 
are in plateaus. 

We present here the elementary landscape decomposition of QAP. In the next 
section we present the formal definition of QAP. Section [3] presents the main 
result and its proof. Finally, we conclude in Section |4] with some conclusions 
and future work. 

2 Quadratic Assignment Problem 

Let P be a set of n facilities and L a set of n locations. For each pair of locations 
i and j, an arbitrary distance is specified nj and for each pair of facilities p and 
q, a flow is specified w pq . The Quadratic Assignment Problem (QAP) consists 
in assigning the facilities of P to the locations in L in such a way that the 
total cost of the assignment is minimized. Each location can only contain one 
facility. For each pair of facilities the cost is computed as the product of the 
weight associated to the facilities and the distance between the locations in 
which the facilities are. The total cost is the sum of all the costs associated to 
each pair of facilities. One solution to this problem is a bijection between P and 
L, that is, x : P —> L such that x is bijective. Without loss of generality we can 
just assume that P = L = {1, 2, . . . , n} and each solution a; is a permutation 
in S n , the set permutations of {1,2,..., n}. The cost function to be minimized 
can be formally defined as: 



n 




(2) 



2 



3 Decomposition of QAP 



In the previous section we defined the search space and the objective function. 
In order to completely define a landscape we need to define the neighborhood. 
The neighborhood N considered here is the swap or 2-exchange neighborhood, 
in which two solutions are neighboring if one can be obtained from the other 
one by a swap (exchange of two elements) in the permutation. In this section 
we prove that the QAP is composed at most by three elementary components 
and we give an expression for them. 

Let us start rewriting (J2j) - In order to analyze the elementary components 
of the fitness functions associated to a problem class, like QAP, it is useful to 
separate in the formal definition of the objective function the information that is 
particular of a given instance (the data of the instance) from the general issues 
that characterize the class of the problem. In the case of QAP, the information 
related to the particular instance is included in the distance matrix (ry) and 
the weight matrix (w pq ). The question now is: how to separate the data of 
the instance in the fitness formulation. There arc different ways to do it, but 
we are interested in linear combinations of functions where the coefficients of 
the functions are associated to the particular instances. The reason for this is 
that any linear combination of elementary functions (with the same character- 
istic constant k) is also an elementary function. With this idea in mind, it is 
not difficult to see that Equation ([2]) can be written using the following linear 
combination: 

n n 

■f^ = EE r » w P9 6 l(i) 6 lu) ( 3 ) 

where we used the Kronecker's delta. The problem-related part of the fitness 
function is, thus, the product ^u\^(jy At this point we can go further and 
deal with a more general objective function. In ([3]) the value of the product 
Tijixipq depends on i, j, p, and q in a particular way, but is not the most general 
one. Using multilinear algebra concepts, the previous product is a four-rank 
tensor that has been computed as a tensor product of two two-rank tensors 
(matrices), which is a special case of four-rank tensor. In the most general 
case we can define a four-rank tensor to replace the product. Let us call i/Jijpq 
the new general four-rank tensor and let us define the parameterized function 
tp(i,j),(p,q)(3i) = ^x(i)^x(j)' Then we can rewrite the fitness function as: 

n 

f = E $Up<lV{i,i),<p,<l) ( 4 ) 

and we can focus our analysis on <{>{i,j),{p,q) ; since any result on it can be extended 
to any linear combination of <p functions, and, thus, to /. Now, the objective 
function of the QAP is just a particular case of our new objective function /, 
in which ipij pq = rijW pq . 
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If i = j and p = q in the functions <f>(i.j),( P .q), then we have f{i,i),(p>p)^ x ^ ~ 
■ For this particular case we have the following 

Lemma 1. Considering the swap neighborhood the function fu^jp^p) is an 
elementary landscape with k = n. 

Proof. In the following, for the sake of clarity we will remove all the parameters 
from the name of the function when there is no confusion. The function ip 
is elementary if and only if there exist two constants a and b such that the 
following expression holds for all the solutions: 

avg{^(y)} = aip(x) + b 

j/£JV(x) 

In order to reduce the expressions we multiply the previous expression by the 
size of the neighborhood, which is d — "^" 2 ~ 1 - > . We then obtain: 

^2 <fi(y) = cip{x) + e (5) 

y£N(x) 

where c = ad and e = bd. Now, we compute the exact expression of J2 y eN(x) fiv) 
for the two different values that tp can take: 

• Case (f(x) = 1 (in this case x{i) = p). From the neighboring solutions 
there are n — 1 with <fi(y) — and the remaining neighbors have a value 
(p(y) — 1. Then we can write: 

J2 <P(j/) = (d-n+l) 

y£N(x) 

• Case ip(x) = (in this case x(i) ^ p). From the neighboring solutions 
there is only one with ip(jj) = 1. The remaining neighbors have a value 
(p(y) = 0. Then we can write: 

E ^y) = l 

y£N(x) 

Now we use Equation (J5]) to obtain the following linear equation system: 

1 1 \ f c \ _ f d-n + 1 
1 J { e ) \ 1 

The solution of the previous system is c = d — n and e = 1; so we have 
a = 1 — n/d and b — l/d. Then, we can write 

&vg{(p(y)} = <p(x) + - = <p(x) + ~ ( - - <p(x)) (6) 

yeN(x) \ dJ d d \n J 

and we conclude that is an elementary landscape with k = n and 

average <P(i,i),(p, P ) = 1/n. □ 
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Before proving the main result of this section we need to introduce a family 
of auxiliary functions that map permutations to R: 



(i,j),(p,q) 



(x) 



a 




= ph x(j) = q 


P 


iix{i) 


= q A x(j) = p 


7 


nx{i) 


= P® x(j) = q 


£ 


if x(i) 


= q © ar(i) = P 


c 


if x(i) 


^p,qA x(j) ^p,q 



(7) 



where 1 < i,j,p,Q < n are integer values with i ^ j and p ^ q and a, /3, 7, £, C € 
R. All the previous values are parameters of the family of functions. We denote 
with © the exclusive-or operator. The previous functions are valuable thanks 
to the following 



Lemma 



a,/3,7,e,C 



is an 



2. Considering the swap neighborhood the function <jj 
elementary landscape in the following cases: 

• a — n — 3, (3 = 1 — n, 7 = —2, e = 0, £ = —1 with k = 2n 

• a = n — 3, P = n — 3, 7 = 0, e = 0, ( = 1 with k = 2(n — 1) 

• a = 2n — 3, P = 1, 7 = n — 2. £ = 0, £ = —1 im'f/i k = n 

Proof. In the following, for the sake of clarity we will remove all the parameters 
from the name of the function when there is no confusion. The function (f> is 
elementary if and only if there exist two constants a and b such that the following 
expression holds for all the solutions: 

avgj>(?/)} = a(j)[x) +b 

y£N(x) 

In order to reduce the expressions we multiply the previous expression by the 
size of the neighborhood, which is d= "( n ~ 1 ) . We then obtain: 



E 

2/£JV(x) 



cf)(y) = aj>{x) + e 



(8) 



where c — ad and e = bd. 

We distinguish five different cases which are symbolically represented in Fig- 
ure [TJ In the figure, each node represents the set of solutions for which one of 
the five branches in ([7]) is true. We label the nodes with the value that <f> takes 
for all the solutions in that node. There exists an arc (i, j) if all the solutions in 
node i have at least one neighboring solution in node j. The label of arc is 
the number of neighbors that any solution in i has in j . Now, we compute the 
exact expression of X^yejv(x) ^(v) f° r the five different values that (j) can take: 

• Case 4>{x) — a. In this case x(i) = p and x(j) = q. From the neighboring 
solutions there is one with cf>(y) = P and 2(n — 2) solutions with <p(y) = 7. 
The remaining neighbors have a value (j>(y) = a. Then we can write: 

4>{y) = P + 2(n-2)j+(d-2n + 3)a 

y£N(x) 
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Figure 1: Transition graph for functions ^i /^fpq)- 

• Case 4>(x) — p. In this case x(i) = q and x(j) — p. From the neighboring 
solutions there is one with <fi(y) = a and 2(n — 2) solutions with 4>{y) = e. 
The remaining neighbors have a value <fr(y) = P- Then we can write: 

<f>(y)=a + 2(n-2)e + {d-2n + 3)P 

y£N(x) 

• Case 4>{x) — 7. In this case x(i) = p or x(j) = q, but not both. From 
the neighboring solutions there is one with 4>{y) = a, two neighbors with 
4>{y) = s, and n — 3 neighbors with 4>{y) = £■ The remaining neighbors 
have a value 4>{y) — 7. Then we can write: 

Hy) = a + 2e+(n-3)C+(d-n) 1 

y£N(x) 

• Case (f>(x) = e. In this case x(i) — q or x(j) = p, but not both. From 
the neighboring solutions there is one with 4>{y) — P, two neighbors with 
4>{y) = 7, and n — 3 neighbors with 4>{y) = (. The remaining neighbors 
have a value 4>{y) = s. Then we can write: 

^2 cj){y) = p + 2 7 + (n - 3)C + (d - n)e 

• Case 4>{x) = (. In this case x(i) 7^ p,q and x{j) ^ p, q. From the 
neighboring solutions there are two with 4>{y) — 7 and two neighbors with 
4>{y) = s. The remaining neighbors have a value <j>(y) = C- Then we can 
write: 

Y, ^(y)=2 7 + 2e+(d-4)C 
Now we use Equation {8J to obtain the following linear equation system: 
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p 

1 

£ 

V c 



i \ 
i 
i 
i 

i / 



/ /3 + 2(n-2) 7 - 
a + 2(n-2)e- 
a + 2e + (n - 
P + 2 7 + (n - 

V 



- (d - 2n - 
-\d-2n- 
3)C + (d- 
3)C+(d- 



3)a \ 

-3)/3 
n)7 



2 7 + 2e + (d - 4)C J 

The previous system has hve equations and two variables, c and e, so it could 
be unsolvable. However, the system can be solved for some value combinations of 
a, P, 7, e, C- In particular, the system can be solved for the value combinations 
mentioned in the statement, that is: 

1. a = n - 3, P = 1 - n, 7 = -2, e = 0, C = -1 

2. a = n-3, /3 = n-3, 7 = 0, e = 0, C = l 

3. a = 2n - 3, P = 1, 7 = n - 2, e = 0, C = -1 

This does not mean that these are the only combinations of parameter values 
for which the system can be solved. They are just three combinations of special 
interest for the goal of this section. It should be noticed here that the linear 
system does not depend on the values of i, j, p, and q. Thus, the solutions to 
the system are also independent of the values of the mentioned parameters. 

Let us study the values of a, b 7 c, e for the first parameter combination, that 
is, a = n — 3, P = 1 — n, 7 = —2, e = 0, and ( = —1. The solution of the linear 
system is c = 5 - ) and e = — 2n, and, thus: a = (1 — 2n/d) and b = —2n/d. 
In order to simplify the notation, let us define ^ .,) ( pg ) = ^^j) (pq)~ ' '~ • 
Then, we can write 



av g{^(ij),( M )(y)} 

yeN(x) 



2n 
~d 



. , 2n / 



-jj n (i,j),(p,q)( X ) - — 



(9) 



and we conclude that fi^ f ^ ^ ^ is an elementary landscape with k = 2n and 
average 0|.. )( ) =-1. 

Let us now focus on the second parameter combination, that is, a = P = 
n — 3,7 = e = 0, and £ = 1. The solution of the linear system is c = (HzlEnzf) 
and e = 2(n — 3), and, thus: a — 1 — 2(n — l)/d and 6 = 2(n — 3)/d. Introducing 



the notation ^f itj )^ p . q ) 



,n-3 : n-3, 0,0,1 



we can write 



av g{^,j),( P , 9 )(2/)} = 

yeN(x) 

= (1 - 2(n - (x) + 2(n - 3)/d = 



- ^,J),(P,9)( X ) + 



2(n-l) 



fi (ij),(p,9)( X ) 



(10) 
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and we conclude that Cth / q % is an elementary landscape with k = 2(n — 1) 



and %UP, q ) = ( n ~ 3 )/(" - 

Finally, let us analyze the third parameter combination, that is, a = 2n — 3, 

/3 = 1, 7 = n — 2, e = and £ = — 1. The solution of the linear system is 



c = "^" 2 3 ^ and e = n, and, thus: a 

i ,. n a ,2n-3,l,Ti— 2,0,-1 

notation -n , \ — <Pr -\ \ 



= 1 — n/d and b 
we can write 



n/d. Introducing the 



avg{^ 3 



(i,j),(p,q)(v)} 
yeN{x) 



n 



= ft 



(11) 



and we conclude that ft 3 4 / ^ is an elementary landscape with k — n and 
O 3 

"(i,j),(p,9) 



1. 



□ 



Now, we arc in conditions of presenting the main result of the paper, which 
is the following 

Theorem 1. For the swap neighborhood defined above, the function f defined 
in is the sum of at most three elementary landscapes with constants k\ = 2n, 
k'2 = 2(n — 1). and k% = n. 



functions ^\ iij)>{p>q y and n ti,j),(p,q) when i ^ 3 and P ^ 1 as ( we omit 



Proof. The functions <P(i,j),(v,q) defined above can be written using the auxiliary 
functions ft}. w >, ft?. ., 

(^,3),{p,q) , (1,3), 

the subindices for clarity) 

1 -ft 2 + , - ^ n 3 (12) 



i 



2n 



1 



2(n - 2) n(n - 2) 

This can be easily appreciated with the help of Table [TJ 



Condition 


u l 


ri 2 




£ 


In 


2(n-2) 


n(n-2) 


x(i) = p A x(j) = q 
x(i) — q A x(j) = p 
x(i) = p © x(j) = q 


n—'A 


n—o 


2n — 


i 





2n 

1-n 


2(n-2) 
n— 3 


n(n-2) 


2n 

-2 


2(n-2) 



n(n-2) 
n-2 


2n 


n(n-2) 


x(i) = q © x(j) = p 














x(i) ^p,q/\x(j) ^p,q 


-1 


1 


-1 





In 


2(n-2) 


n(n-2) 



Table 1: Elementary components and their sum. 



Since the ft family of functions are elementary, the ip family of functions are 
a sum of three elementary components, namely: j^ft 1 , 2 (n-2) and n(n-2) 
This decomposition of 93 allows us to write the fitness function / as a decompo- 
sition of elementary landscapes in the following way: 
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/ = ''Pijpq V(i,i),(p,g) + V'iipp ^(iiOiCPiP) = 

i, j, p, ij = 1 i,P=l 

p + q 

n /q1 Q2 o3 \ 

_\^,/, f "(»j).(p,g) | "(».j).(p,g) , "(j,j),(P,g) \ 
-2^Wm I 2n + 2(n-2) n(n-2) ) 

i + 3 
rt 

"i" ^ ' ^iipp V(m)>(p>P) 
z,p— 1 

The elementary components of / are: 
n SI 1 

/d= E^^f^ 1 (13) 

i,j,P, 9=1 
i # 3 
P # 9 

™ SI 2 

i # J 
P # 9 

n £^3 n 

/c3 = 1 P 2 ) + E ^"PP ^(i.i).(p.P) ( 15 ) 
p, q = 1 i.p=l 
« ^ J 
P # 9 



where the functions / c i, / C 2, and / C 3 are elementary with constants k\ = 2n, 
k2 = 2(n — 1) and k% = n, respectively, because they are a linear combination 
of elementary functions. Thus, / can be written in a compact form as / = 

fcl+fc2+fc3- ' □ 

In the statement of the theorem we say that the number of elementary 
components is three at most. That is, this number cannot be larger than three, 
but it could be lower. It is possible that for some particular instances the number 
of elementary landscapes could be reduced (we will see later that this happens 
for the TSP). 

Unlike elementary landscapes, when the fitness function is the sum of several 
elementary components with different characteristic constants, the average value 
in the neighborhood does not linearly depend on the fitness function. Since, f c \ , 
f C 2, and / C 3 are elementary components, the Grover's wave equation ([1} can be 
applied to them. And we can compute the average value in the neighborhood 
in the following way: 
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avg{/(y)} = avg{/ c i(y)} + avg{/ c2 (y)} + avg{/ c3 ( 2 /)} 

y£N(x) y£N(x) y£N(x) y€N{x) 

= fcl(x) + f c2 (x) + f c3 (x) + ^ (/ cl - f cl ( X )) + 

+ ^ {fc2 ~ fc2(x)) + | (/c3 - U(X)) = 

= f{x) + -i- (/ c i - / c i(z)) + - (/ c2 - / c2 (z)) 
n — 1 n 

+ — — r (/c3 - /ca(a:)) 
n — 1 

4 Conclusions and Future Work 

We have proven that QAP can be decomposed as a sum of three elementary 
landscapes when the swap neighborhood is used (Theorem [T]). We have pre- 
sented the exact expressions for the three components in Equations (flU)) , (TT^j) . 
and (fT5|) . The elementary components of QAP allow one to compute the av- 
erage value of the objective function in the neighborhood of a given solution x 
using the evaluation of the three components in x. 
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